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CHAPTER  1 

INTRODUCTION 

The  idea  of  constructing  a  very  large  optical  aperture  hy 
combining  several  smaller,  independent  optical  systems  is  a  relatively 
new  one  to  the  optics  community.  Although  this  synthesis  is  standard 
for  radar  systems,  the  precision  required  to  phase  optical  wavelength 
posed  a  seemingly  unattainable  requirement.  Recent  demands  for  verv 
large  telescopes  to  advance  astronomical  frontiers  and  to  perform 
defense  missions  have  inspired  innovative  soiutinns  to  this  prohl“m. 
Within  the  past  few  years,  the  idea  of  aperture  synthesis  has 
established  itself  as  the  most  promising  technique  for  constmcino 
mammoth  optical  systems. 

Optical  phasing  of  synthetic  aperture  systems  require  that  each 
separate  aperture  have  the  same  optical  character) sties  with  regard  tn 
far  field  pattern  and  optical  path  length. n  In  order  to  phase 
multiple  telescopes  together  one  must  ensure  that  the  wavefronts,  as 
they  are  leaving  their  respective  telescopes,  are  absolutely  phased 
with  respect  to  one  another.  They  must  also  he  kept  pointing  tn  the 
same  spot  in  the  optical  far  field.  The  Phased  Array  (PHASARl 
Experiment  is  one  such  effort  and  is  an  on-going  experiment  at  the  Air 
Force  Weapons  Laboratory,  Kirtland  AF8,  Mew  Mexico.  »'  figure  1 
illustrates  the  experimental  layout  of  the  PHASAR  system. 
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Figure  1.  Experimental  Layout  of  PHASAR  System. 


^T|y»|ji  ywi^jw'1 1  >«■  p^w  p%j*  i  ■••  f^<«p 


This  experiment  consists  of  three  telescopes  mounted  within  a  single 

housing.  Each  telescope  has  a  dedicated  set  of  optical  path 

difference  adjusters.  These  adjusters  steer  a  flat  mirror  so  that 

each  adjuster  is  capable  of  changing  the  optical  path  length,  ti’t 

error,  and  phase  error  of  its  respective  telescope.  These  adjusters 

must  he  coordinated  in  a  common  reference  system  so  t*»at  continuous 

phasing  can  be  achieved.  A  number  of  electrical,  optica1,  and 

electro-optical  devices  are  being  utiTized  to  address  these  con+rni 

issues.  Peripheral  control  issues  include  open-loop  control  of  the 

adjusters  to  allow  for  phasing  of  the  telescopes  w^ile  one  or  more 

apertures  are  moved  in  piston,  and  global  determination  of  the  point 

where  there  is  no  opticcal  path  difference  between  apertures. 

Figure  2  shows  the  component  block  control  diagram. 

* 

This  project  is  being  performed  within  prescribed  limits 
determined  by  the  thesis  committee  and  the  experimenters.  Firstmost, 
the  experimental  optics  are  considered  to  be  ideal  hy  the  models 
considered  here.  These  optics  have  been  previous!y  modeled  by  and  are 
well  documented  in  PHASAR  technical  memorandums,  so  here  the  focus  of 
the  effort  is  toward  the  remainder  of  the  experiment.  Additiona'Hv, 
the  tracker  leg  of  this  experiment  has  been  excluded  from  modeling. 
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Figure  2.  Component  Block  Control  Diagram. 


CHAPTER  2 


MODELING 

Referring  again  to  Figure  2  we  see  that  the  block  diagram  is 
composed  of  two  types  of  control  loops:  Tilt  and  Piston.  This 
chapter  outlines  the  modeling  process  for  these  loops  in  two  major 
subsections. 


2.1  TILT  LOOP 

The  tilt  loop  was  modeled  with  six  components  as  shown  in 
Figure  3.  The  two  inputs  and  two  outputs  model  the  X  and  Y  positions 
of  the' beam  centroid  striking  the  quad  cell.  There  are  three  tilt 
loops  in  the  system,  one  controlling  each  aperture  in  tilt. 
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Figure  3.  Tilt  Loop  Block  Diagram. 


2.1.1  COMPENSATOR 

The  compensator  was  built  by  the  experimenters  in  order  to 
operate  the  tilt  error  control  loops  within  a  desired  bandwidth, 
schematic  of  the  compensator  is  shown  in  Figure  4.  The  companion 


notch  filter  was  inserted  at  the  point  marked  A.  It  is  shown  in 
Figure  5.  The  compensator  can  be  manually  switched  from  a  low  pass 
filter  to  an  integrator.  It  is  normally  operated  as  an  integrator  in 
closed-loop  and  as  a  low  pass  filter  during  open-loop.  The  component 
transfer  functions  were  modeled  after  these  schematics.  '  Figure  fi 
contains  a  Bode  plot  of  this  component  model. 

Integrator  TF  : 

s  +3.2673els  +4.0012e2s'+?.1767e3s~+4.4-We3s 

,  TC  s4+1.5708els3+1.2337e2s2+1.046Be3s+4.a3fl2e^ 

Low  pass  TF  : - 3 - * - n - 

( l+F.305e-3s)  (  s  +  3.2673els  M.0012e2s  +?.1767e3s+4.438?e^ 

2.1.2  AXIS  SEPARATOR 

The  axis  separator  is  used  to  transform  the  signals  which  command 
the  Optical  Path  Difference  Adjustor  (OPDA),  to  induce  x-tilt,  y-tilt 
and  piston  corrections  to  the  apertures.  The  axis  separator  is 
implemented  in  various  places  via  hardware  summers  and  gain  pots.  Thn 
reason  for  the  axis  separator  is  that  misalignment  errors,  the 
imperfect  positioning  of  the  Piezo-Electric  Tranducers  (PZTs)  on  the 
OPDAs  and  the  coupling  between  the  PZTs  represent  errors  that 
necessitate  some  means  of  minimization. 

Misalignment  errors  are  the  results  of  not  having  the  x-y 
coordinate  frame  of  the  quad  cell  oriented  identically  to  the 
corresponding  coordinate  frame  of  the  OPOA.  Figure  7  illustrates  this 
point. 
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Figure  5.  Notch  Filter  Schematic. 
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Figure  6.  Notch  Filter  Component  Model  Bode  Plot. 
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Figure  7.  Coordinate  Misalignment  Axis. 


The  axis  separator  can  he  heuristically  seen  as  compensating  for  the 
misalignment  angle  by  modifying  the  movement  commands  sent  to  the 
PZTs.  A  simple  direction  cosine  matrix  will  perform  this  function. 
With  this  misalignment  angle  nominally  very  small,  the  sine  terms  are 
at  least  two  orders  of  magnitude  less  than  the  cosine  terms. 
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The  positioning  of  the  P7Ts  on  the  OPDAs  is  the  major  reason  for  the 
axis  separator.  It  must  transform  commands  in  terms  of  x  and  y  into 
individual  commands  to  the  three  PZTs,  dl,  d?  and  d3,  as  shown  in 
Figure  8. 
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Figure  8. 

PZT  Placement  on  Actuator . 

A  simple  translation  matrix  perfonxs  this  function.  Here  the  pain  Is 
set  at  2.5. 

r dl“l  f  n.o  -l.o]  r COS  a  sin 

d2  *  2.5  j  -1.0  0.5 C0S 

Ld3j  L  1-0  0.5J 

2.1.3  HIGH  VOLTAGE  AMPLIFIER 

The  high  voltage  amplifier  (HVA)  has  a  single  stage  gain  of  200 

and  was  custom  built  for  this  experiment  by  Los  Alamos  National 
Laboratory.  They  are  used  to  amplify  the  signal  sent  to  the  PZTs  in 
the  OPOAs.  Figure  9  represents  a  Bode  plot  of  a  typical  HVA.  The 
main  features  in  this  figure  are  the  unifom  rolloff  of  ?Mb  per 
decade  and  zero  phase  lag.  Taking  the  average  of  the  plots  from  the 


n 


12  different  HVAs,  the  corner  frequency  was  calculated  to  be  1.458 
KHz.  The  single  pole  model  with  the  given  corner  frequency  proved 
more  accurate  than  a  second  order  model  that  was  heavily  damped.  So 
the  transfer  function  for  this  component  was  chosen  as: 

200 

1  ♦  1.0916e-4s 


2.1.4  ACTUATOR 

The  optical  path  difference  adjusters  (OPDA),  or  actuators,  were 

built  for  this  experiment  by  the  Rocketdyne  Division  of  Rockwell 

International.  The  actuator  is  a  dynamic  flat  mirror  that  can  be 

positioned  with  great  precision.  ‘Three  PZTs  are  housed  in  a  flexture 

1  1 

that  forms  a  mount  for  the  mirror.  When  a  PZT  is  actuated,  it 
physically  pushes  against  the  back  of  the  mirror,  moving  it  in  tilt  or 
translating  it  in  piston.  If  all  three  PZTs  are  actuated  equally,  the 
mirror  translates  in  piston  with  little  tilt  movement.  Since  the  pZT 
has  a  dynamic  range  of  only  five  microns,  motorized  micrometers  were 
stacked  on  the  rear  of  each  PZT  to  provide  coarse  adjustment  in  open 
loop  operation.  A  Bode  plot  was  drawn  for  the  actuator  and  is  shown 
in  Figure  10.  The  dominant  features  of  this  plot  are  the  134  db  per 
decade  rolloff  and  the  resonant  behavior  in  the  vicinity  of  2.7  KHz. 
The  OPOA  was  visibly  vibrating  near  this  frequency.  The  gain 
magnitude  must  be  discounted  because  it  was  adjusted  so  as  to  fit 
within  the  physical  limits  of  the  plotting  equipment.  The  model 
excludes  the  resonant  behavior  because  it  lies  well  outside  the 
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bandwidth  of  the  present  experiment.  After  many  discarded  attemots,  a 
seventh-order  model  was  developed  to  represent  this  behavior.  A  Bode 
plot  of  Its  transfer  function  Is  shown  In  Figure  11.  The  gain  of  this 
transfer  function  corresponds  to  the  plus  or  minus  BOH  volts  from  the 
HVA  translating  to  a  range  of  plus  or  minus  2.5  microns.  ?o  the 
transfer  function  for  this  component  is: 


5.0e-d 
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The  x  or  y  displacement  of  the  quad  cell  is  a  function  of  the  angle  e 
and  the  attributes  of  the  optical  setup.  Here,  M  *  S'/S  «  40.0,  and 
F  a  0.3  meters.  Then  the  displacement  ad  *  FeM  or  ad  *  12®.  An 
expression  for  e  was  developed  from  the  geometries  of  the  actuator 
shown  in  Figure  13. 


Figure  13.  OPDA  Geometry. 

To  determine  ax  a  right  triangle  is  formed  by  the  distance  a,  which 
was  determined  to  be  0.01386  meters,  and  the  commanded  di sp1 acement. 
This  triangle  is  shown  in  Figure  14. 
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Figure  14.  Tilt  Optics  OPDA  Triangle. 


Here,  9  *  arctan(  d3/a) ,  and  hence  ax  a  12  arctan(  d3/a) .  Similarly,  we 
have  &y  »  12  arctan( -dl/c) .  Since  the  arctan  function  in  this  rangp 
can  be  linearized^  to  arctan(k)  *  k,  our  component  transfer 
functions  are: 
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2.1.6  QUADRANT  CELL  AND  UDT  431 

The  quad  cell  is  used  to  determine  the  displacement  of  the  beam 
centroid  striking  its  surface  from  the  coordinate  origin.  The 
detector  used  in  this  experiment  is 'a  UDT  Spot./AD,  position  sensing 
photo  detector  quadrant  cell.  See  Figure  IS. 
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Figure  IF,  puaq  Cell  Cross  Section. 
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The  quad  cell  sensor  has  a  bandwidth  of  1  MHz  and  is  therefore  far 
outside  the  observed  closed-loop  system  bandwidth  of  around  1  KHz. 

The  UOT  431  X-Y  position  monitor  is  a  processing  unit  that  has  the 
individual  signals  A,  B,  C  and  0  as  inputs  the  X  and  Y  di splacements 
along  with  other  information  not  presently  used  in  this  experiment  as 
outputs.  The  position  equations  for  this  detector  use  A,  B,  C  and  0 
to  represent  electrical  terminal  signal  outputs.  The  equations  are: 

„  „.  ,  .  '  (A  ♦  D)  -  (B  ♦  C) 

X  displacement:  x  =  ~b " ♦  "PT" 

v  ..  .  *  (A  ♦  B)  -  (C  ♦  D) 

Y  displacement:  x  =  —  -g- 

From  existing  test  data  on  the  UDT  431,  it  was  found  that  the  X  and  Y 
displacement  outputs  have  a  bandwidth  of  about  1  KHz.  This  is  caused 
by  stages  A21  and  A25  driving  both  the  position  output  as  well  as  an 
analog  meter.  This  limits  bandwidth  due  to  the  high  open-loop  output 
impedence  of  the  operational  amplifier  driving  a  788  ohm  meter  as  well 
as  the  output  load.  This  low  value  of  load  resistance  reduces  the 
open-loop  operational  amplifier  gain  and  as  a  result  limits  the 
bandwidth  potential  of  the  gain  of  10  (open-loop)  output  driver  (see 
Figure  16).  This  test  data  also  verified  that  the  analog  divide 
module  output  was  slew  rate  limited  to  1.6  KHz.  The  gain  would  be  a 
function  of  the  signal  processing  if  it  were  not  for  the  analog  divide 
stage  normalizing  the  output.  So  then  using  the  X  and  Y  displacement 
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equations,  gain  is  a  function  of  rli spl acement,  spot  size,  intensity 
and  detector  area.  As  previously  stated,  the  detector  is  2. 67mm 
square  and  the  diameter  of  the  beam  striking  the  cell  Is  2.337mm.  Tt 
has  constant  intensity.  Then  the  gain  can  be  strictly  described  by 
the  displacement  equations  as  long  as  the  beam  does  not  lap  over  the 
edges  of  the  detector  area.  Although  the  tilt  loop  has  the  dynamic 
range  to  drive  the  beam  off  the  detector  area,  this  event  is 
considered  catastrophic  and  the  result  of  an  external  disturbance. 

The  gain  of  the  UDT  431  corresponding  to  the  entire  beam  in  the  right 
half  plane  (RHP)  is  10  volts,  and  -10  volts  for  the  entire  beam  in  the 
left  half  plane  (IHP).  What  was  needed  to  model  this  behavior  was  an 
equation  relating  centroid  displacement,  and  beam  radius,  to  outDut 
qain.  If  the  beam  area  is  completely  in  the  RHP,  then  the  output,  gain 
can  be  modeled  as: 

Gai n  *  ( 2A ‘ - irr? ) 

it  r 

So  now  an  equation  is  needed  for  A'.  >'sing  the  approximation  *rom 
Figure  17  and  proceeding  classically.  A'  can  be  represented  as: 


A’  ,  it  r  - 
r-h  s  ad 
a- 

A  *  it  r  - 

2 

s  Tr 

A’  =  irr  ‘  - 


[•“  -  l(r-h)] 


C-4  -  lad] 


-  2iad|V  r2-ad?l 
r*Cos’lM  -  2  Adi  V7V 

r  3 
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Figure  17.  Quad  Cell  Gain  Geometry 


Now  combining  this  with  the  gain  equation  gives: 


„  .  lOad  2,2  -l’adL 

Gain  =  — - - r( 2*r  -2r  cos  —■  J  + 

wr  Ad  r 


4  [adl 


2  -.2 
fr  -Ad 


However,  given  that  the  entire  beam  remains  on  the  detector  area,  the 
maximum  ax  and  Ay  can  be  determined  from  comparing  the  beam  diameter 
and  detector  size.  It  becomes  apparent  that  a  value  greater  than 
0.16^nm  for  ax  and  Ay  will  displace  the  beam  centroid  such  that  a 
portion  of  the  beam  laps  off  the  detector  area,  therefore,  ax  and  av 
have  a  nominal  range  of  plus  or  minus  0.169mm.  Examining  the  behavior 
of  this  gain  equation  over  this  small  nominal  range,  it  was  found  that 
the  following  rectangular  approximation  was  still  more  accurate.  Let: 


2 

A'  =  2rad  ♦  ±j- 

2 

Gain  =  [2(2rAd  *  -  * r?] 


then, 

^Oa  d 
irr 


So  with  the  maximum  ax  or  Ay  *  0.169mm  then  the  error  is  less  than 
4.4e-9  square  meters  or  the  relative  error  is  less  than  0.174%.*' 
This  figure  is  quite  acceptable  and  this  approximation  is  nicely 
linear.  The  overall  transfer  function  for  the  quad  cell  is  then: 


where  K  is  a  loop  sensitivity  factor  of  0.1. 


2.2  PISTON  LOOP 

The  piston  loop  was  modeled  in  six  components  and  is  shown  in 
Figure  18.  The  deviation  with  respect  to  the  reference  position  is 
modeled  by  a  single  input  and  single  output  model.  There  are  also 
three  piston  loops  in  the  model,  two  control ing  a  single  aperture  in 
piston  and  a  third  serving  as  a  reference.  The  component  models  from 
the  HVA  and  actuator  are  identical  to  those  used  in  the  tilt  loop. 


Figure  18.  Piston  Loop  Block  Diagram. 

2.2.1  COMPENSATOR 

The  piston  loop  compensator  differs  from  the  tilt  loop 
compensator  only  by  a  gain  factor  of  0.2.  This  was  necessary  to  step 
down  the  gain  from  the  phase  management  data  processor  to  plus  or 
minus  one  volt  for  the  axis  separator. 

2.2.2  AXIS  SEPARATOR 

In  the  piston  loop  the  axis  separator  model  is  very  similar  in 
function  to  the  tilt  loop  model,  except  that  here  it  must  onlv 
interpret  the  movement  conwands  for  one  incut.  *  simole  translation 
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matrix  will  translate  a  piston  movement  command  into  individual  PZT 
displacement  commands. 
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2.2.3  PISTON  OPTICS 

As  previously  stated,  the  optical  transmission  characteristics 
are  outside  the  scope  of  this  effort,  but  could  be  combined  with  this 
model  if  a  joint  study  is  desired.  Here  the  behavior  modeled  was  the 
translation  of  a  given  displacement  from  the  actuator  into  an  optical 
path  difference  (OPD).  This  transfer  function  is  a  function  of  the 
weighted  average  of  the  six  different  wavelengths  lasing  in  the 
experiment;  here,  x  =  0.5  microns.  It  can  he  seen  in  Figure  19- that 
this  transfer  function  is  periodic  every  2x.  The  nominal  operating 
range,  save  a  gross  external  disturbance,  is  between  -0.F  microns  onn 
and  *0.5  microns  OPD.  Therefore  the  transfer  function  i$  simply  -].Q 
within  the  nominal  operating  range. 

2.2.4  PHASE  MANAGEMENT  OATA  PROCESSOR 

The  phase  management  data  processor  (PMPP)  was  built  for  this 
experiment  by  the  BDM  Corporation. 12  It  is  composed  of  three  major 
subsystems:  A  Fairchild  CCD  1300  Line  Scan  Camera  (LSC),  and  camera 
control  circuitry;  analog-to-digital  converter  (ADC);  comouter 
processor.  The  aperture's  beam  is  interfered  with  the  reference  bean 
and  the  resulting  interferogram  light  pattern  is  focused  on  the  charge 
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coupled  device  (CCD).  This  interferogram  is  sampled  every 
millisecond,  passed  through  the  ADC,  and  is  represented  as  a  linear 
array  of  256  pixels  with  nine  bits  per  pixel.  Figure  20  shows  an 
interferogram  that  has  been  imaged  onto  the  LSC.  This  plot  was  made 
by  outputting  a  single  frame  of  digital  data  from  the  PMDP.  The 
central  portion  of  the  digitized  pattern  is  used  for  estimating  the 
OPD  between  the  two  beams. ^  The  algorithm  for  a  linear  measure  of 
0P0  is: 

OPD  * 

where  R  is  the  right  side-lobe  and  L  is  the  left  side-lobe.  This 
algorithm  compares  the  integrated  intensity  of  the  two  adjacent  side 
lobes  and  normalizes  the  quantity  by  their  sum.  Simulated  results  of 
this  algorithm  are  plotted  in  Figure  22.  Under  computer  control  *he 
actuator  was  moved  a  given  step  size  while  taking  data  from  the  LSC 
and  with  the  PMDP  calculating  the  OPD  values.  Figure  21  was  produced 
from  this  data.  The  widening  of  the  plot  toward  the  right  side  was 
caused  hy  the  nonlinearity  of  the  old  HVAs.  Recause  of  the 
ambiguity,  it  is  not  possible  with  this  algorithm  to  know  when  the  two 
path  lengths  are  exactly  equal,  corresponding  to  zero  fringe.  A 
global  algorithm  is  used  to  determine  zero  fringe  while  the  system  is 
in  open-loop.  This  linear  algorithm  is  then  used  for  closed-loop 
control.  The  modeling  of  this  behavior  was  divided  into  two  parts. 

The  optical  translation  of  displacement  into  OPP  was  used  to  develop  a 
transfer  function  for  the  optics  in  the  piston  loop  (see 
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Section  2.2.3).  OPD  is  then  the  input  to  the  PMDP,  with  a  gain  in 
volts  to  correct  to  zero  OPD  as  the  output.  Figure  23  Illustrates 
this  function  of  PMDP  gain.  To  model  the  one  millisecond  delay,  an 
e‘Ts  term  proceeds  the  transfer  function  for  the  PMDP.  A  zero-order 
hold  is  included  to  form  a  piecewise  continuous  output  signal,  and 
therefore,  a  (l-e"^s)/s  term  is  cascaded  onto  the  transfer 
function.  In  both  instances,  T  =  0.001  seconds. 

2.3  TILT-PISTON  COUPLING 

Referring  again  to  Figure  2,  there  are  two  coupling  mechanisms 
that  potentially  interweave  the  piston  and  tilt  loop  of  each 
aperture.  These  are  descriptively  called  piston  corrupting  tilt  (PCT) 
and  tilt  corrupting  piston  (TCP).^  PCT  has  been  observed  in  this 
experiment  but  had  little  significance  as  it  was  30db  down  from  the 
uncoupled  tilt  signal.  TCP  is  potentially  significant  because  it 
could  corrupt  the  piston  reference  and  thus  couple  together  the  three 
apertures.  However,  during  the  operation  of  this  experiment,  the 
magnitude  of  this  coupling  was  never  significant  enough  to  allow 
measurement.  So  coupling  in  this  model  will  be  variable  in  magnitude, 
but  whose  significance  will  be  along  the  magnitude  previously  stated. 

During  the  alignment  phase  of  the  experiment,  the  axis  separator 
is  adjusted  via  the  gain  pots  and  summers  such  that  the  cross  coupling 
signals  are  at  least  30db  down  from  the  desired  signal.  This  can  be 
heuristically  seen  as  aligning  the  misalignment  angle  so  that  the 
cross  coupling  between  individual  apertures  is  insignificant. 


31 


CHAPTER  3 


REALIZATION 

The  next  step  in  the  modeling  scheme  is  to  develop  dynamical 
equations  from  the  transfer-function  matrices.  For  every  realizable 
transfer-function  matrix  there  is  an  unlimited  number  of 
dynamical -equation  real i zations.  Therefore  a  major  consideration  in 
the  realization  process  is  to  find  a  "good"  realization.  Here  a 
realization  that  was  both  controllable  and  observable  was  pursued.^ 

3.1  TILT  LOOP 

Combining  the  component  models,  in  cascade,  from  the  previous 
chapter  a  forward-loop  transf er-function  (FLTF)  matrix  was  obtained. 
The  matrix  was  unit  scaled  to  Krad/sec  to  keep  the  coefficients 
manageable.  This  is  shown  in  Figure  24.  Referring  back  to  Figure  2, 
a  closed-loop  transfer-function  (CLTF)  matrix  with  unity  feedback  was 
computed  and  is  shown  in  Figure  25.  The  CLTF  matrix  was  transformed 
into  an  observable  canonical  form.  The  modal  matrix  corresponding  to 
this  system  was  badly  conditioned  and  precluded  many  realization 
approaches.  It  was  determined  that  a  high-order,  reducible 
realization  would  suffice.  This  realization  form  is  a  decoupled 
composite  dynamical  equation^  and  is  illustrated  in  Figure  26.  The 
realization  can  be  easily  accomplished  by  inspection,  but  a  penalty  is 
incurred.  This  realization  is  not  internally  coupled  and  the 


33 


FLTF  Tilt  Loop. 


rare 


an(s)  a12(s) 
9n(s)  g'J2(s) 

a21(s)  a22{s^ 

g^TTT  g^TsT 


with  ajj(s)  defined  in  Figure  24. 


gn(s)  -  s15  +  8.6314*la14  +  3.5499*3s13  +  9.340le4a12  +  l^SSeda11  +  2.5561e7a10 
+  2.9089e8s9  +  2.6462e9s8  +  1.9317elQs7  +  1.1243ells6  +  5.1245ella5 
+  1.7715el2s4  +4.3971el2s3  +  7.144lel2a2  +  7.8914el2s  +  1.1938el3 

gl2(s)  -  a15  +  8.6314eU14  +  3.5499e3a13  +  9.3401e4a12  +  1.7685C6811  +  2.5561e7s1(3 
+  2.9089e8s9  +  2.6462e9sa  +  1.9317elOs7  +  1.1243eUs6  +  5.1245eUs5 
+  l.7702el2a4  +  4.3765el2s3  +  6.9824el2s2  +  6.5I93eL2a  +  6.1186*12 

gn(s)  •  a15  +  8.6314ela14  +  3.5499 *3a13  +  9.3401e4a12  +  1. 768566a11  +  2.5561e7s10 
+■  2.9089e8s9  +  2.6462e9s8  +  1.9317el0a7  +  1.1243ells6  +  5.1245ells5 
+  1.7688«12a4  +  4.3544el2a3  r  6.8094el2s2  +  5.0522el2s  -  1.0356ell 

g22(a)  -  s15  +  8.6314ela14  +  3.5499*3a13  +  9.340ie4a12  +  1.7685e6s11  +  2.556Le7a10 
+  2.9089«8a9  +  2.6462e9a8  +  1.93l7elOs7  +  1.1243*IU6  +  5.1245ells5 
+  1.771UL2*4  +  4.39H*12s3  +•  7.10i2eL2a2  +  7.5273el2s  +  1.0393*13 

Figure  25.  CLTF  Tilt  Loop. 
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Figure  26a.  Tilt  Loop  Realization. 


Figure  26b.  Tilt  Loop  Realization. 
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dimension  is  undesirably  large.  Model  reduction  will  be  an  important 

follow-on  issue.  The  "CONTROL .lab"  computer  package,  made  available 

through  the  Laboratory  for  Computer-aided  Design  of  Systems  and 

Robotics,  was  used  extensively  during  this  process  and  proved  itself 
1 ft 

invaluable. 


3.2  PISTON  LOOP 

Combining  the  component  models,  in  cascade,  from  the  previous 
chapter  a  FLTF  was  obtained.  It  was  again  unit  scaled  to  Krad/sec  to 
remain  consistent  with  the  tilt  loop,  and  to  keep  the  coefficients 
manageable.  This  is  shown  in  Figure  27.  Referring  back  to  Figure  2, 
a  CLTF  with  unity  feedback  was  computed  and  is  shown  in  Figure  28.  A 
realization  was  obtained  by  solving  the  14th  order  differential 
equation  for  the  highest  derivative  term.  It  was  then  transformed 
into  a  nested  sequence  of  integrations.  From  this  form  a  simulation 
diagram  was  drawn  and  the  state  differential  equations  obtained  (see 
Figure  29).  This  procedure  leads  to  an  irreducible,  observable 
canonical  form  real ization. *  This  is  shown  in  Figure  30. 

3.3  STABILITY 

The  experimenters  were  particularly  interested  in  the  stability 
of  the  model,  since  this  was  a  system  attribute  that  they  were  exposed 
to  in  their  daily  experimental  work.  The  tilt  loop  was  chosen  for 
closer  examination  because  of  greater  experimenter  familiarity  and  the 
relative  ease  of  inspection.  A  pole-plot  of  the  FLTF  is  illustrated 
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Figure  28.  CL TF  Piston  Loop. 
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Figure  29.  Piston  Loop  Realization  Procedure. 
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Figure  30.  Piston  Loop  Realization. 


in  Figure  31;  the  poles  are  tabulated  in  Table  1.  The  pole  at  the 
origin  represents  the  integrator  contained  within  the  compensator.  A 
pole-plot  of  the  corresponding  CLTF  matrix  is  shown  in  Figure  32,  with 
poles  tabulated  in  Table  2.  Both  of  these  plots  were  broadly  accepted 
by  the  experimenters  as  representative  of  known  system  behavior.  The 
marginally  RHP  poles  in  the  CLTF  pole-plot  correspond  to  the  constant 
coefficients  in  terms  g^,  g,,^,  and  gj£  from  the  CLTF  matrix. 

The  experimenters  claim  that  this  instability  does  indeed  exist  prior 
to  open-loop  adjustments  made  at  the  onset  of  an  experimental  run. 

3.4  LARGE-SCALE  SYSTEM 

The  combined,  three-aperture  model  that  corresponds  to  Figure  2 
consists  of  three  tilt  loops  and  two  piston  loops.  A  realization  of 
this  208th  order  system  is  shown  in  Figure  33  using  the  tilt  and 
piston  loop  matrix  partitions  defined  in  Figures  26  and  30 
respectively.  It  is  apparent  that  this  system  is  coupled  only  through 
the  functions,  f  and  e.  As  was  stated  earlier,  although  the  exact 
mathematical  characteri sties  of  these  functions  are  not  known,  their 
general  dynamics  are  bounded  above  by  a  constant  M,  where  M  is  30db 
down  in  magnitude  from  the  pure  signal.  The  magnitude  of  e  is  still 
smaller,  and  has  never  been  significant  enough  to  facilitate 
measurement. 
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Table  2. 


Tilt  CLTF  Matrix  Poles. 


A1 1  Poles  A12  Poles 


Real  Part 

Imam' narv  Part 

Real  Part 

Imaginary  °art 

♦ 

0.2486 

- 

1.5939 

- 

0.1451 

- 

1.3085 

•f 

0.2486 

♦ 

1.5939 

_ 

0.1481 

♦ 

1.3086 

- 

1.9016 

- 

6.0063 

_ 

2.0866 

- 

8.4pqa 

- 

1.9016 

6.0063 

- 

2.0866 

♦ 

8.488° 

- 

2.0050 

- 

8.6034 

- 

2.3657 

- 

6.8272 

- 

2.0050 

+ 

8.6034 

- 

2.3657 

+- 

6.8270 

- 

3.7082 

- 

8.9987 

_ 

3. 52°5 

_ 

0.8611 

- 

3.7082 

♦ 

8.9987 

- 

3.52^5 

♦ 

8.8611 

- 

7.9799 

- 

0.6418 

_ 

7.8557 

_ 

0.8°63 

- 

7.9799 

+ 

0.6418 

. 

7.8557 

0 . 806° 

- 

8.2251 

- 

6.0505 

- 

7.°84? 

- 

r.4?l t 

- 

8.2251 

♦ 

6.0505 

- 

7.984? 

♦ 

5.4214 

- 

11.0270 

- 

3.8851 

- 

10.6374 

- 

7.°681 

- 

11.0270 

♦ 

3.8851 

- 

10.6374 

♦ 

3.9681 

- 

17.1175 

- 

0.0000 

- 

17.1057 

- 

0.0000 

A21  Poles  A22  Poles 


Real  Part 

Inar 

jinary  Part 

Real  Part 

Imag-' 

nary  Pa 

♦  0.2000 

♦ 

0.0000 

0.1548 

- 

1.5336 

-  2.1454 

- 

8.2588 

+ 

0.1648 

♦ 

1.5336 

-  2.1454 

♦ 

8.2588 

• 

2.0040 

_ 

5.9663 

-  3.1923 

- 

8.8154 

- 

2.0040 

•f 

5.3*63 

-  7.1923 

♦ 

8.8154 

2.0260 

- 

8.5711 

-  4.0693 

• 

5.2807 

- 

2.0260 

♦ 

8.5711 

-  4.0693 

♦ 

5.2807 

- 

3.6599 

- 

9.°634 

-  4.3571 

- 

2.2317 

- 

3. 65°9 

♦ 

8.9634 

-  4.3571 

♦ 

2.2317 

- 

7.8596 

- 

n.oooo 

-  4 , 94  94 

♦ 

0.0000 

- 

8.1316 

- 

5,0740 

-  8.2201 

- 

2.3058 

- 

8.1316 

♦ 

c . 9340 

-  8.2201 

♦ 

2.3058 

- 

8.1936 

♦ 

O.oon° 

-  10.1617 

. 

4.5045 

• 

10.317? 

- 

7,0844 

-  10.1617 

♦ 

4.5045 

- 

10.317? 

9, OQ4  4 

-  17.0327 

+ 

0.0000 

_ 

17.1130 

0.0000 

Figure  33.  Large-Scale  System  Realization. 


CHAPTER  4 


CONCLUSIONS 


A  large-scale  realization  of  this  experimental  system  has  heen 
obtained.  This  model  incorporates  the  potential  for  comolete  couolinq 
between  tilt  and  piston  loops  of  the  same  aperture,  although  all  data 
indicates  that  the  system  is  weakly  coupled.  As  with  most  projects  of 
this  type,  as  many  questions  have  been  raised  as  conclusions  drawn. 
Further  examination  of  the  coupling  could  be  particularly  interesting, 
and  applied  to  other  multiple  aperture  systems  wi.th  similar 
attributes.  Individual  component  models  could  be  extended,  or  used  in 
their  present  form,  for  use  in  other  system  modeling  efforts.  The 
combination  of  this  model  with  the  separate  system  optical 
characteri zation  could  be  most  illuminating.  This  modpl  represents  a 
valuable  investigative  tool  for  the  experimenters  to  explore  various 
control  issues  such  as  increased  system  bandwidth,  and  n<»ar-0Dtimum 
control.  The  development  of  more  robust  control  algorithms,  that 
would  continue  to  provide  asymptotic  stability  and  regulation  under 
perturbation  of  plant  parameters,  could  also  be  explored.  A  priori 
information  could  also  be  obtained  regarding  a  potential  hardware 
component  upgrade.  The  model  could  be  expanded  to  include  the  DC 
drift  tracker  leg  to  more  completely  duplicate  the  entire  PWA^AR 
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experiment,  or  extended  and  generalized  for  use  as  a  design  tool  for 
multiple  aperture  systems. 

This  work  also  represents  a  point  of  departure  for  additional 
issues  that  could  warrant  further  study.  A  full-scale  computer 
simulation  would  decisively  verify  the  modeling  of  this  system  and 
would  provide  a  computational  aid  for  the  experimenters.  The 
simulation  should  be  accomplished  with  inputs  set  to  zero,  and  driven 
only  by  a  disturbance  injected  into  the  quad  cell  and  PMDP.  This 
disturbance  rejection  simulation  will  most  closely  match  experimental 
conditions.  Model  reduction  would  be  of  great  value  for  a  system  this 
large  in  dimension.  In  particular,  the  reduction  of  the  tilt  loop  to 
a  more  manageable  order  would  significantly  reduce  the  complete 
model's  order.  There  are  a  number  of  model  reduction  schemes  that 
could  be  applied  to  this  system;  the  balanced  method  appears 
particularly  attractive.14  This  approach  uses  concepts  of 
controllability  and  observability  to  determine  a  reduced-order  model. 
The  full-order  model  has  many  unobservable  or  uncontrollable  modes, 
and  is  therefore  an  appropriate  candidate  for  this  method. 

Near-optimum  control  techniques  could  be  applied,1^  and  the  system 
decomposed  via  interaction  prediction  to  find  a  near-optimal 
hierarchical  control  regime.  These  efforts  would  further  strengthen 
thi s  model ' s  utility. 
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APPENDIX 


GLOSSARY  OF  ABBREVIATIONS 


AOC  -  Analog  to  Digital  Converter  -  A  signal  converter  that  is  a 
subsystem  of  the  PMDP. 

CCD  -  Charge  Coupled  Device  -  An  1024  element  photo  sensor  used  as  an 
image  detector  in  the  LSC  subsystem  of  the  PMDP. 

CLTF  -  Closed-Loop  Transfer  Function. 


FLTF  -  Forward-Loop  Transfer  Function. 

HVA  -  High  Voltage  Amplifier  -  A  key  component  in  the  PHASAR 
experiment. 

LHP  -  Left  Half  Plane  -  Usually  used  in  reference  to  the  s-plane 
during  stability  considerations. 

LSC  -  Line  Scan  Camera  -  An  electronic  camera  system  that  produces  an 
analog  video  output  waveform  representing  the  spatial 
distribution  of  the  time-integrated  brightness  across  one 
scanned  line  of  an  optical  image.  A  subsystem  of  the  PMDP. 

OPD  -  Optical  Path  Difference  -  The  difference  in  path  length  between 
two  apertures. 

OPDA  -  Optical  Path  Difference  Adjuster  -  A  key  component  in  the 
PHASAR  experiment. 


PCT  -  Piston  Corrupting  Tilt  -  The  phenomenon  of  commanded  piston 
displacements  inducing  uncommanded  tilt  displacements. 

PHASAR  -  PHASed  ARray  experiment  -  The  experiment  of  interest  located 
at  the  Air  Force  Weapons  Laboratory. 


PMDP  -  Phase  Management  Data  Processor  -  A  key  component  in  the  piston 
loop  of  the  PHASAR  experiment. 

PZT  -  PieZo-electric  Transducer  -  The  key  dynamical  element  in  the 
OPDA. 


RHP  -  Right  Half  Plane  -  Usually  used  in  reference  to  the  s-plane 
during  stability  considerations. 

TCP  -  Tilt  Corrupting  Pi ston  -  The  phenomenon  of  commanded  tilt 
di splacements  inducing  uncommanded  piston  displacements. 


UDT 


United  Detector  Technology  -  A  position  sensing  photodetector 
manufacturer.  Their  equipment  is  a  key  component  in  the  tilt 
loop  of  the  PHASAR  experiment. 
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